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Abstract: Knowledge of i^e-Fe/Pb differential cross sections for i/^ energy below several tens of MeV scale 
is believed to be crucial in understanding Supernova physics. In a segmented detector at Spallation Neutrino 
Source, energy reconstructed from the electron range measurement is strongly affected because of both 
multiple scattering and electromagnetic showers occurring along the electron passage in target materials. In 
order to estimate the effect, a simulation study has been performed with a cube block model assuming a perfect 
tracking precision. The distortion of energy spectrum is observed to be proportional to the atomic number of 
target material. Feasibility of unfolding the distorted Ve energy spectrum is studied for both Fe and Pb cases. 
Evaluation of statistical accuracy attainable is therefore provided for a segmented detector. 
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1 Introduction 



Precise knowledge of the differential cross sections 
for quasi-elastic neutrino-nucleus scattering at a scale 
up to several tens of MeV is of importance in under- 
standing the explosion dynamics of supernova, since 
the neutrinos play a decisive role in the development 
between the prompt shock and the delayed shock 
Among the cross sections neutrino interactions with 
^H, C, O, Fe and Pb are most relevant. In this low 
energy scale, the cross section varies rapidly, leading 
to large uncertainty in predicting the core collapse 
supernova [1, 0] ■ Although some theoretical calcula- 
tions 01 of neutrino-nucleus cross sections have been 
performed in recent years, large uncertainties still ex- 
ist and need to be clarified from experiments. On the 
experiment side, the differential cross sections with C 
and O have been measured by KARMEN and LSND 
experiments Q , but those with Fe and Pb are still 
not available. In this paper, we present a simulation 
study on cross section measurement of neutrino inter- 
action with Fe and Pb. The motivation is to give an 



estimate of the achievable limit in the experiment of 
neutrino cross section measurement with a segmented 
detector at Spallation Neutron Source(SNS). 

At SNS, protons are driven to hit a target, produc- 
ing a large amount of neutrons and by-product pions. 
The negative pions are mostly captured in the target, 
while the positive pions, after having been stopped 
in the target within 0.1 ns or so, undergo the char- 
acteristic successive decay scheme 7r+ "'^> /i+ + 1^^ 
and followed by "°> e+ + J^e + S'm ■ The resul- 

tant j/j's can have a flux up to 2 x 10^ /cm^/s at a 
distance of 20 meters away from the spallation target 
for the SNS at the Oak Ridge National Laboratory, 
assuming the SNS power to be 1 MW 01 • Since the 
time structure of SNS beam provides a unique advan- 
tage for neutrino studies 0, E3: Ell; backgrounds are 
strongly suppressed. In the SNS target, the pions are 
stopped and decay to muons. The subsequent muon 
decay produces equal intensity of and j/,,, with a 
maximum energy 52.8 MeV. The energy spectra of 
the three kinds of neutrinos are shown in FiglT] 

Such an energy range overlaps extremely well with 
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that of supernova neutrinos. This characteristics 
makes it valuable to measure cross sections of interac- 
tions between neutrinos and various target materials 
at SNS. Fig. [5] shows the sketch of neutrino detector 
at SNS. A J/— SNS project [tI] was proposed to build a 
neutrino facility at the SNS at Oak Ridge National 
Laboratory. However, there is no report on how pre- 
cise the measurement can be achieved. 
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Fig. 1. Neutrino spectrum from muon and pion 
decay-at-rest at a spallation neutron source. 
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Fig. 2. 
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Sketch of neutrino detector at SNS. 



Electron neutrinos interacting with target mate- 
rial at the SNS undergo both neutral- and charged- 
current interactions, while muon neutrinos can only 
undergo neutral-current interactions due to the lim- 
itation of kinematics. The neutral-current events 
are suppressed by at least two orders of magnitude 
in comparison with the charged-current events and 
can only be measured by the decay of excited nu- 
cleus. The charged-current events can be detected via 
the outgoing electron in the neutrino-nucleus reaction 
i/e+ zXjv — ^+iYjv-i+e~ (i^e +ri ^ p + e"). Since 
the neutrino energy is much lower than the nucleus, 
considering the charged-current quasi-elastic(CCQE) 
process and ignoring fermi motion of nucleus, the neu- 
trino energy E^,^ can be determined by the induced 
electron energy i?^. 



(1) 



where A is the mass difference between the final state 
nucleus Y and the initial state nucleus X. There- 



fore, a precise energy reconstruction of the i/^, in- 
duced electrons is essential to determine the differ- 
ential cross sections. The segmented detector is pro- 
posed to measure the trajectories of induced elec- 
trons in the target. However, the energy lose of 
the electron is dominated by ionization at low en- 
ergy, while above the critical energy E^, it will be 
dominated by bremsstrahlung. The electromagnetic 
shower caused by bremsstrahlung will lead to large 
uncertainty on the measurement of electron track 
length for energy above Ec- The critical energy can 
be approximated |81, 191] by E, = 800 MeV/(Z 1.2), 
where Z is the atomic number of the target nucleus. 
Therefore the critical energy of Fe is E^ = 29.4 MeV, 
and Ec = 9.6 MeV for Pb. This makes it impossible to 
obtain the electron energy by its track length event by 
event. In this paper, we will show a statistical appli- 
cation to unfold the distorted energy spectrum, and 
leading to a reasonable resolution of the differential 
cross section. 

The outline of this paper is as followed. In section 
2, the Geant4 simulation of neutrino induced elec- 
trons are simulated in Fe and Pb. The simulation 
results and analysis are described in section 3, and a 
brief conclusion is made in section 4. 

2 Simulation procedure 

To simplify the issue, we use a cube block tar- 
get material with a fiducial volume of 1 x 1 x 1 m^, 
in stead of the segmented detector with target, tube, 
wire, etc. The tracking resolution is assumed to be 
perfect. Geant4 (version 9.0p01) package 1^ is ap- 
plied in this simulation study. 

The cross section of neutrino-nucleus scattering is 
rather small, we therefore directly simulate the be- 
haviors of the induced electron in the target. The 
input energy spectrum of electron is the Michel dis- 
tribution shaped by the theoretically predicted cross 
section given by R. Lazauskas and C. Volpe The 
maximum energy of induced electron, -Emax, is ap- 
proximately 48.2 MeV on ^^Fe, and the threshold en- 
ergy of the interaction is E'thrcs — 5.1 MeV. For ^°*Pb, 
E^,^ ~ 49.9 MeV, and E.t,,,, ~ 3.4 MeV. 

Daughter products e.g., ^^Co, ^°*Bi are ignored, 
since the life times of the daughters are long compar- 
ing the search window (~ 10 micro seconds). 

The material of the detector is selected to be iron 
or lead with natural abundances, but not pure ^^Fe 
or ^"*Pb. Therefore, the measured cross section is an 
average of mixture isotopes of iron or lead. This af- 
fects iron little since the natural abundance of ^''Fe is 
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91.754%. For lead, the natural abundances of ^''^Pb, 
207pb and 206pb are 52.4%, 22.1% and 24.1% respec- 
tively, the effect of the mixture should be considered 
seriously in later studies. 

For the interaction of electrons in nuclear mate- 
rial, electromagnetic process is the most important 
one. The processes of multiple scattering, low energy 
ionization and low energy bremsstrahlung are added 
for electrons in the simulation. 

To perform the unfolding procedure l3"15|, two 
kinds of spectra of electrons are generated: 1) Michel 
distribution used to obtain the detector response ma- 
trix; 2) Michel distribution weighted with the cross 
section of corresponding z^e-nucleus interaction. The 
second one is considered as the true distribution to 
be reconstructed. For each event, the electron's en- 
ergy (input) E and track length L in the target are 
recorded. 

The number of events attainable in one nominal 
year (assuming 10^ seconds) is estimated by consid- 
ering the flux of electron neutrino, the size of detector 
and the average cross section of I'e-nucleus scattering. 
The expected maximum cross section at 52.8 MeV for 
Fe is cTpcmax — 1.2 X 10~^^ cm^, and it gives that the 
average cross section over energy is about O.Scrpcmax- 



For Pb, cTpe,, 



1.3 X 10 cm^, and the average 



cross section over energy is about 0.4apb,max- Assum- 
ing the effective volume of target is 1 m^, and the 
flux of 1/5's is 2 X 10^ cm~^s~^, the number of events 
in one nominal year is estimated to be 6.1 x 10^ for 
Fe target, and 3.4 x 10'' for Pb target. 



3 Unfolding method 

In the energy measurement, due to detector reso- 
lution, efficiency and sometimes intrinsic physical rea- 
son, the true value in bin i, will usually migrate into 
measured energy bin j. In other words, the measured 
distribution is usually distorted by the detector re- 
sponse, or the detector response is folded to the true 
distribution, leading to the measured distribution. In 
many cases, the true distribution could be restored 
by applying correction factors to the measured dis- 
tribution. This method fails, however, if the relation 
between true value and measured ones is highly non- 
linear where the behavior of detector response is poor. 

Suppose the true distribution is Xtmc, and the 
measured one is l^bs- If effects of efhciency and back- 
ground are ignored, they satisfy 

Yobs — R' X^^^c, (2) 

where R is the response matrix or the migration ma- 



trix which describes the effects of detector response. 
The basic idea of unfolding is that the effects of detec- 
tor response could be removed or unfolded by inverse 
of the response matrix, i.e., the estimator of the true 
distribution X is 



^true — R ' Y^l 



(3) 



However, although this is the unbias estimator with 
smallest variance, the physics is totally swept out by 
unphysical oscillation. The regularization is needed 
to remove the unphysical oscillation, which is actu- 
ally a compromise between bias and unphysical oscil- 
lation, controlled by the choice of regularization pa- 
rameter. There are several different schemes of reg- 
ularisation (l3l - [l5| . for example, Tiknonov reulariza- 
ton, regularization based on entropy, Bayesian reg- 
ularization, regularization based on Singular Value 
Decomposition(SVD), etc. The basic point of all is 
that the true distribution is believed to be smooth 
but not oscillating rapidly. In this study, we used the 
TSVDUnfold class in ROOT packagelSj. TSVDUn- 
fold is a class in ROOT for unfolding based on SVD 
technique, it is part of RooUnfold, a ROOT unfolding 
framework, developed by Tim Adye, et al [l^ . 

In the SNS neutrino experiment, to measure the 
differential cross section of neutrino-nucleus scatter- 
ing, the true value is the energy of incident electron, 
which is principally unknown and supposed to be 
binned in histogram E, the measured value is the elec- 
tron track length in the target, binned in histogram 
L. The response matrix R can be obtained by Monte 
Carlo simulation as long as we have full knowledge on 
the detector. In discrete expression. 



^j=^^ii-E.i {j = l,...,nt), 



(4) 



where Ei is number of events in bin i of E, Lj is 
number of events in bin j, n„ and Uf, are numbers of 
bins of E and L respectively. The physical meaning 
of Rji is the probability of an event in bin i of true 
distribution to migrate into bin j of measured dis- 
tribution. Therefore, the response matrix is actually 
the normalized two-dimensional scattering plot of L 
and E. 

4 Simulation results and analysis 

In the simulation, the electron track length L is 
directly read from Geant4 with no uncertainty. The 
response matrix of detector with the iron target is 
shown in Fig. [3ja), and that of detector with the 
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lead target is shown in Fig. |3l[b), corresponding to 
200,000 electrons with Michel distributed energy. If 
the matrix were diagonal, one can easily obtain the 
energy of incident electron through track length event 
by event. However, these two plots show significant 
non-diagonal elements of the response matrices, and 
hence it is impossible to obtain directly the incident 



energy through track length. Especially, it is obvious 
that the track length in iron target is smeared widely 
with energy above 30 MeV, while for lead target the 
smearing begins with energy above 10 MeV or so. 
This is the consequence of electromagnetic shower. 
Both simulations are exactly consistent with the ex- 
pected critical energies of Fe and Pb. 
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Fig. 3. Scattering plot of electron energy E versus track length L for target of Fe (a) and Pb (b). Significant 
non-diagonal elements of response matrix exist due to electromagnetic shower above critical energy Ec- The 
energy spectrum of incident electrons is Michel distribution. 



To extract the physical part out of the unphysi- 
cal oscillation, the choice of regularization parameter 
is crucial, which reflects a trade-off between bias and 
oscillation. For the unfolding based on the SVD tech- 
nique, the choice of regularization parameter should 
be tuned for given distribution, number of bins and 
sample size. As described in Ref. 117|, the regulariza- 
tion parameter kterm should be chosen by the plot 
of log|c?|i versus i, where integer i is from 1 to the 
number of bins and the ith component of the vec- 
tor d is the coefficient in the decomposition of the 
measured (and re-scaled) histogram. Only the first 
few terms, say fc, of the decomposition should be sig- 
nificant for a reasonably smooth measured distribu- 
tion, while the others correspond to contribution of 
quickly oscillating basis vectors, and should be com- 
patible with zero. Therefore, one should see two sep- 
arate patterns on the plot of log|d|i: for small i, \di\ 
is significantly greater than 1 and falls gradually to a 
standard gaussian distribution for large i. One should 
choose the regularization parameter kterm to be the 
critical value i = k, after which d^'s are insignificant. 

Fig. m^a) shows the plot of log|c?i| versus i with 
5,000 events accumulated in Fe target, approximate 
data of one nominal year. The parameter kterm — 5 
is chosen from this plot. The unfolded energy spec- 
triun of the v. induced electron is shown in the upper 



plot of Fig. Elja), while the lower plot shows the dis- 
crepancy between the unfolded spectrum and the true 
one, where the latter is Michel distribution shaped 
by theoretical prediction of differential cross section. 
With the unfolded energy spectrum of the induced 
electron, one can obtain the spectrum of the 's that 
interact with the Fe target using Eq. ([T]). Dividing 
this spectrum by Michel distribution gives the differ- 
ential cross section of i^e-nucleus scattering, as shown 
in the upper plot of Fig. [B^a). The lower plot in 
Fig. m^a) shows the relative discrepancy between the 
measured cross section and the true one. One can see 
that the discrepancies and errors are huge for energy 
below 20 MeV, mainly because the cross section of 
low energy is too small and few events are collected 
in this range. The discrepancy of the last bin is also 
large, around 30%. For other bins, the discrepancies 
are controlled within 20%. We checked this result by 
comparison with five other samples under the same 
condition. It shows that for energy below 20 MeV the 
measured cross sections ffuctuate significantly, there- 
fore totally unreliable. While the discrepancies be- 
tween 20 and 50 MeV are well controlled within 20%. 
Increasing the sample size to be 10000 slightly im- 
proves the behavior between 20 and 50 MeV, but 
helps less (if not nothing) to the region below 20 MeV. 



No. X YANG Zhen-wei et al: Simulation Study on neutrino nucleus cross section measurement in Segmented Detector 



5 




2 4 6 8 10 12 14 16 18 20 



' 2 4 6 8 10 12 14 16 18 20 



(a) (b) 

Fig. 4. Plots of log|di| versus i for target of Fe (a) and Pb (b). The number of bins is 20 for both targets, 
while the number of measured events is 5,000 for Fe and 30,000 for Pb. 
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Fig. 5. Comparison between the unfolded and true energy spectra of the induced electrons. The target is 
Fe for (a) and Pb for (b). The number of bins is 20 for both targets, while the number of measured events 
is 5,000 for Fe and 30,000 for Pb. 
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Fig. 6. Comparison between the measured and true differential cross sections of z/e-Fe (a) and i^e-Pb (b) 
scattering. The number of bins is 20 for both targets, while the number of measured events is 5,000 for Fe 
and 30,000 for Pb. a(E ^52.8MeV) is 1.2 x 10"^''cm^ for Fe and 1.3 x 10"^**cm^ for Pb. 



For Pb target, the sample size is 30,000. Fig.lDJb) 
shows the plot of log \di\ versus i, according to which, 
the regularization parameter kterm = 4 is chosen for 
Pb. The unfolded energy spectrum of the induced 



electron is shown in the upper plot of Fig. El^b), and 
the lower plot shows the discrepancy between the 
unfolded and true spectrum, where the true one is 
Michel distribution shaped by theoretically predicted 
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cross section of J^e-Pb interaction. The resulting dif- 
ferential cross section of z^e-Pb scattering is shown in 
the upper plot of Fig. Eljb), and the lower plot shows 
the relative discrepancy between the measured differ- 
ential cross section and the true one. Similar to the 
case of Fe target, one can see that the discrepancies 
and errors are large for energy below 22 MeV due to 
lack of events. The discrepancy of the last bin is also 
large, around 30%. For other bins, energy between 
25 and 50 MeV, the discrepancies are well controlled 
within 20%. This result is compared with five other 
samples using the same regularization parameter. It 
shows that for energy below 22 MeV the measured 
cross sections fluctuate significantly, therefore totally 
unreliable. While the discrepancies between 22 and 
50 MeV are well controlled below 20%. 

5 Conclusions 

With Monte Carlo simulation, we studied the 
achievable accuracy for the J/e-nucleus scattering cross 
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